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SUMMARY 2 $£>3^ 

Correlations for variable property heat-transfer and friction pressure 
drop data were obtained for forced convection of hydrogen and nitrogen through 
electrically heated tungsten wire mesh. These correlations represent the data 
of six different helically coiled wire meshes for the following range of condi- 
tions : 

(1) Mesh porosity of 64 to 72.2 percent 

(2) Wire diameter of 0.020 to 0.035 inch 

(3) Surface temperature of 1400° to 5200° R 

( 4 ) Outlet gas temperature of 600° to 2400° R 

(5) Mass velocity for hydrogen of 0.4 to 3.1 lb/(sec)(sq. ft) and for 

nitrogen of 4.5 to 10.2 lb/(sec)(sq. ft) 

(6) Heat flux of 0.5 to 8.3 Btu/( sec )(sq. in.) 

(7) Pressure level of 1 atmosphere 

The effect of flow bypass, resulting from a mesh heater not filling the flow 
passage, was investigated on a 0.030-inch-diameter wire mesh for a bypass area 
of 25 percent. 

INTRODUCTION 

The Lewis Research Center is conducting research on a tungsten-water- 
moderated, hydrogen propelled, nuclear rocket concept. A part of the experi- 
mental phase of this program is to evaluate various types of fuel elements and 
supporting structures at simulated operating conditions. 




The results of this report were used to design the heating elements of 
a high- temperature hydrogen preheater for hot flow testing. The prerequisites 
for these elements were that they must have: 

(1) Capability of being electrically heated to surface temperatures of 

5500° R 

(2) Electrical resistance to match an available high-voltage power supply 

(3) Sufficient surface area to transfer the generated heat to the flowing 

gas 

(4) Sufficient flow area to minimize gas pressure drop 

Commercially available tungsten mesh made of interwound helical coils of 
tungsten wire were considered applicable as the heating elements for this pre- 
heater, but their heat-transfer and pressure drop characteristics were not 
known. 

A literature survey revealed that there is a limited amount of experimen- 
tal heat-transfer and pressure drop data available for forced convective flow 
through porous wire mesh (referred to by other authors as porous media) and the 
majority of these data are for constant property conditions. Reference 1, 
which summarizes the results of the Stanford-Office of Naval Research program 
on compact heat-transfer surfaces, provides most of the existing data on in- 
dividual wire mesh elements of different wire diameters and porosities. By 
using a transient test technique and constant property conditions, this program 
determined heat-transfer correlations for each mesh. Reference 2 revised these 
data and obtained a general correlation for all mesh. Similarly the isothermal 
pressure drop data for each mesh correlated individually with the Fanning type 
equation, but a general correlation for all meshes was not obtainable. Refer- 
ence 3 provides the only available variable property heat-transfer and pressure 
drop data. The data were obtained for steady-state flow of air through an 
electrically heated tube bank at surface temperatures up to 1100° R. 

A considerable amount of pressure drop data are available for fluid flow 
through packed beds (ref. 4). A comparison between wire mesh and packed beds 
indicates the pressure drop characteristics of both are dependent on the same 
basic parameters, that is, mass velocity, fluid properties, and geometrical 
factors. The correlations reported for packed beds are used as a basis for 
correlating the pressure drop data of this report. 

This report presents experimental variable property pressure drop and 
heat-transfer correlations for helically coiled wire mesh at average surface 
temperatures to 5200° R. Six different meshes were tested. The wire diameter 
varied from 0.020 to 0.035 inch with mesh porosities between 0.640 to 0.722. 

The mass velocity through these mesh was varied from 0.4 to 3.1 pounds per 
second per square foot for hydrogen and from 4.5 to 10.2 pounds per second 
per square foot for nitrogen. The heat flux ranged from 0.5 to 8.3 Btu per 
second per square inch. 
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APPARATUS AND PROCEDURE 

A schematic diagram of the flow system, test section, power supply, in- 
strumentation, and corresponding components associated with each, as used in 
this investigation, is shown in figure 1. 


Flow System 

As may be seen in figure 1, hydrogen or nitrogen was supplied to the flow 
system from a tube trailer at a maximum pressure of 2400 pounds per square 
inch. From the trailer the gas then flowed through a preset pressure reducer 
valve, a remotely operated control valve, a choked flow nozzle and an on-off 
valve that supplied gas directly to the test section. The gas flow was metered 
by means of the choked flow nozzle that assured a constant mass flow through 
the test section. From the test section the heated gas flowed through a two- 
baffle molybdenum mixing can and into a gas to water concentric tube heat ex- 
changer where it was cooled below 1000° R before being exhausted into the 
atmosphere . 

For safety purposes, the entire system was purged with nitrogen before 
hydrogen entered the system, and the controls were set for fail safe operation 
so if a predetermined safety permissive stops the hydrogen flow, nitrogen 
would automatically purge the system. In such a case, the electrical test 
power would also be automatically shut down. 
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Figure 2. - Tungsten wire mesh. 
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TABLE I. - MESH GEOMETRIC PARAMETERS 


Mesh number 


Wire diameter, d, in. 0.020 0.025 0.030 0.035 0.020 0.025 

Mandrel diameter, D, in. 0.060 0.060 0.080 0.080 0.060 0.060 

Equivalent diameter, D e , ft 0.00425 0.00372 0.00635 0.00642 0.00425 0.00404 

Mesh size (length X width), 3X1 3X1 3X1 3X1 3X3 3X3 

in. X in. 

Number of parallel coils, N 23 23 17 16 69 69 

Coil pitch, p, in. /turn 0.0666 0.0814 0.111 0.125 0.0666 0.077 

Porosity, e 0.722 0.640 0.717 0.685 0.722 0.660 





Figure 3. - Exploded view of test section assembly. 


Test Section 

Experiments were performed on a test section consisting of mesh formed by 
interwound helical coils of tungsten wire as shown in figure 2 . Each end of 
the coils is sandwiched between two tungsten plates, approximately 0.060 inch 
thick, and the coils are heliarc welded at the ends of the plates to provide 
positive mechanical and electrical connections. 

Four different 3- by 1-inch meshes were tested, each having a different 
wire diameter and porosity. Also tested were two 3- by 3-inch meshes, one of 
which has the same parameters as one of the 3- by 1-inch meshes. Table I lists 
the geometrical parameters associated with each mesh. The smaller 3- by 1-inch 
meshes were used for most of the data because they permitted investigation of 
higher ranges of mass velocities and heat fluxes and approximated the geometry 
of the final heater design. 

Figure 3 shows an exploded view of the test section assembly, which con- 
sists of a 3-foot entrance transition section to straighten the gas flow prior 
to entry into the test section, a boron nitride housing designed to hold and 
electrically insulate the mesh bus connections and to minimize bypass of the gas 
around the mesh, a water- jacketed stainless-steel outer housing, and a 
molybdenum can to mix the gas prior to measuring its temperature. A rubber 
0-ring near the cold end of the boron nitride housing prevented leakage of the 
gas between the boron nitride and the stainless-steel outer support housing. 

No provision was made for expansion of the mesh. During initial heating, 
bowing occurred in the flow direction, and the mesh retained a permanent set 
after cooling. This set could not be changed even by reversing the mesh and 
heating to 5000° R at high flow rates. 


Power Supply 

A single phase 60-cycle 500-kilovolt-ampere saturable-reactor controlled 
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power supply was used to electrically heat the tungsten mesh. Output voltage 
was varied from 4.7 to 50 volts with a maximum current rating of 10 000 amperes. 
With bus losses, however, the maximum power to the test element was limited to 
225 kilowatts . 


Instrumentation 

The location of the instrumentation is shown in figure 1. The voltage 
across and the current through the test section were measured. The test 
section voltage was taken directly across the mesh to eliminate any error 
caused by a voltage drop between the power supply and the test section. A true 
root-mean-square voltmeter was used to measure the test voltage because of the 
wave form produced by the saturable-reactor controlled power supply. Current 
was read on a precision ammeter through a 4000:1 step down current transformer. 
Inlet pressure to the test section was measured with a calibrated 0- to 100- 
pound-per-square-inch Bourdon tube gage. The pressure drop across the mesh 
was continuously recorded with a ±l-pound-per-square-inch temperature- 
compensated strain-gage bridge differential pressure transducer. Inlet temper- 
ature was measured with a type K thermocouple (designation ref. 5), and the 
exit temperature was measured with a platinum/platinum- 13-percent-rhodium 
thermocouple. The exit thermocouple was placed in a baffled molybdenum mixing 
can to give a true mixed bulk gas temperature. Pressure and temperature 
measurements at the inlet to the choked flow nozzle were made. The mass flow 
rate was set by adjusting the nozzle inlet pressure. 


METHOD OF CALCULATION 
Geometrical Factors 

The mesh heating elements were made of interwound helical tungsten coils. 
These mesh can be completely specified by five parameters: wire diameter d, 

mandrel diameter D, number of parallel coils N, length of mesh b, and helical 
coil pitch p. (All symbols are defined in the appendix.) 

The geometrical parameters and corresponding equations associated with the 
calculation of the heat transfer, pressure drop, and mesh surface temperature 
for the data of this report are as follows: 

The length of a single helical coil S is given by the equation 

S = | ^(D + d) 2 + p 2 (1) 

while the thickness of the mesh in the direction of flow is given by 

L = D + 2d (2) 
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and the total heat-transfer surface area for W number of coils is 


A s = ndSN 


The equivalent diameter for porous media is normally defined as 


Eh 


4(void volume) 
A a 


4LA. 


fl 


( 4 ) 


* where the average flow area is given as 


A fi ~ ^ft 


(5) 


The porosity is defined as the average porosity for the entire mesh volume by 

( 6 ) 


Volume of mesh - Volume of tungsten in mesh _ 


Volume of mesh 


Average Heat-Transfer Coefficients 

The average heat-transfer coefficient was computed from the experimental 
data by the relation 


W c p,b( T 2 


T i ) 


h = , q r = ••Q ) ^ DV v t 1 

V T s ’ T b ) A s (T s - T^I 


(?) 


where 


T b “ 


T x + T 2 


( 8 ) 


The average surface temperature T s of the mesh was determined from the 
relation between temperature and resistivity of tungsten as given in refer- 
ence 6. The resistivity was calculated by 


^ IS 


(9) 


where 


A e = 




( 10 ) 


Average screen surface temperatures from the resistivity temperature relation 
calculated by equation ( 9 ) were verified by means of an optical pyrometer that 
was sighted through a port located in the inlet transition section of the test 
section housing. Optical pyrometer measurements indicated that the entire mesh 
operated within ±500° R of the same temperature over its entire frontal area. 
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Average Friction Pressure Drop 

Static pressure taps were located in the boron nitride tunnel close to 
the mesh to minimize any pressure drop due to friction on the tunnel walls. 

The total static pressure drop across the mesh AP^ is made up of three 
components, a combined expansion and contraction pressure drop, a momentum 
pressure drop, and a friction pressure drop. 

Neither the expansion nor the contraction losses are experimentally 
separable from the friction pressure drop as in the case of tube-type flow. 
They are therefore included as part of the friction pressure drop. 

The relation between these components is expressed by the following equa- 
tion: 


AP t 


g£ / i_ 

« v 2 


— ) + AP 
p l/ 


( 11 ) 


where G^ /g[(l/p 2 ) - (1/p )] is the momentum pressure drop and AP is the 

friction pressure drop (including expansion and contraction effects). The mean 
density for the friction pressure drop calculations as used in equation (ll) 
is evaluated from the static pressures and total temperatures of the gas as 
follows : 


P 


m 


1 F 1 + P 2 
R T x + T 2 


( 12 ) 


The difference between the total and static temperatures of the gas was never 
more than l/2 percent. Total gas temperatures were therefore used in these 
calculations . 


RESULTS AND DISCUSSION 

The data obtained from this investigation are presented in table II. 
Differences between the total electrical heat input and the heat transferred 
to the gas were due to radiation and conduction losses to the water-cooled 
buses and jacketed housing. 


Heat Transfer 

Flow normal to a wire mesh heat-transfer surface is somewhat analogous to 
flow past a single cylinder in an infinite fluid because the flowing fluid 
forms a laminar boundary layer on the front portion of the cylinders . Proceed- 
ing around a cylinder, the flow accelerates and then decelerates, which causes 
separation of the boundary layer from the surface producing a turbulent wake 
behind the cylinder. The heat transferred from the upstream portion of the 
cylinder where the laminar boundary layer exists can be calculated from the 
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TABLE II. - TABULATION OF HEATED AND ISOTHERMAL DATA 


Hun 

Mesh 

number 

Voltage, 

V 

Current, 

I, 

amp 

Average 

surface 

tempera- 

ture, 

T s > 

°R 

Gas 

Inlet 

gas 

tempera- 

ture, 

T i> 

°R 

Outlet 

gas 

tempera- 

ture, 

? 2 > 

°R 

Gas 

flow 

rate, 

w, 

Ib/sec 

Inlet 

pressure, 

, p i’ 

Ib/sq In. gage 

Total 

static 

pressure 

drop, 

AP t , 

lb/sq in. 

Heated 

1 

1 

28.5 

1280 

2283 

Hydrogen 

520 

985 

0.0199 

2.40 

0.159 

2 



39.1 

1340 

2875 





1185 

.0199 

2.60 

.219 

3 



50.1 

1565 

3110 





1185 

.0297 

5.20 

.395 

4 



39.2 

1520 

2585 





1025 

.0297 

5.00 

.325 

5 



29.5 

1445 

2120 





885 

.0297 

4.70 

.265 

6 



30.6 

1560 

2048 





825 

.0410 

7.20 

.194 

7 



43.6 

1650 

2638 





965 

.0410 

7.80 

.454 

8 



51.5 

1690 

2980 





1070 

.0410 

8.10 

.504 

9 



49.5 

1640 

2950 





1100 

.0349 

6.60 

.435 

10 



29.0 

1500 

2023 





845 

.0349 

5.90 

.305 

11 

2 

15.0 

1805 

1618 

Hydrogen 

520 

870 

0.0199 

2.60 

0.259 

12 



29.7 

1960 

2680 





1235 



3.00 

.409 

13 



22.9 

1880 

2225 





1055 



2.90 

.339 

14 



37.7 

2080 

3152 





1455 



3.30 

.489 

15 



43.8 

2100 

3500 





1630 



3.50 

.539 

16 



49.8 

2160 

3800 





1815 



3.60 

.589 

17 



48.4 

2360 

3450 





1470 

.0297 

6.10 

.835 

18 



39.3 

2280 

2980 





1260 



5.70 

.715 

19 



28.6 

2155 

2400 





1040 



5.20 

.555 

20 



19.4 

2080 

1800 





860 



4.80 

.445 

21 



20.7 

2220 

1780 





810 

.0410 

7.50 

.634 

22 



31.8 

2360 

2435 





985 



8.20 

.804 

23 



41.0 

2470 

2900 





1130 


. 

8.80 

.964 

24 



50.9 

2560 

3345 





1310 


f 

9.40 

1.124 

25 

3 

21.9 

1800 

2460 

Hydrogen 

520 

975 

0.0199 

2.70 

0.169 

26 



31.0 

1900 

3153 





1200 



2.90 

.229 

27 



41.8 

2000 

3870 





1470 



3.10 

.269 

28 



41.0 

2165 

3570 





1235 

.0297 

5.40 

.400 

29 



45.8 

2205 

3850 





1320 



5.60 

.425 

•30 



31.1 

2140 

2950 





1040 



5.10 

.325 

31 



19.7 

1920 

2120 





815 



4.60 

.225 

32 



31.2 

2205 

2790 





925 

.0410 

7.80 

.434 

33 



40.0 

2300 

3330 





1060 



8.20 

.514 

34 



45.5 

2360 

3610 





1150 



8.50 

.564 

35 



17.2 

1960 

1865 





750 

.0349 

5.60 

.245 

36 



28.5 

2110 

2695 





925 



6.20 

.355 

37 



39.5 

2220 

3395 





1110 


. 

6.70 

.455 

38 



45.0 

2280 

3700 





1200 


f 

7.00 

.495 

39 



20.0 

1780 

2300 





925 

.0199 

2.60 

.149 

40 



41.2 

2000 

3820 





1455 

.0199 

3.10 

.269 

41 

4 

23.6 

2170 

2900 

Hydrogen 

520 

1160 

0.0199 

2.80 

0.239 

42 



39.0 

2420 

4020 





1655 



3.30 

.369 

43 



44.3 

2520 

4338 





1845 


T 

3.50 

.409 

44 



42.4 

2680 

3963 





1445 

.0297 

5.90 

.575 

45 



36.2 

2570 

3585 





1290 



5.60 

.515 

46 



23.1 

2390 

2638 





965 


T 

5.00 

.375 

47 



23.0 

2550 

2490 





885 

.0411 

7.60 

.494 

48 



35.0 

2720 

3345 





1110 



8.40 

.654 

49 



44.3 

2860 

3900 





1280 


f 

9.00 

.764 

50 



23.2 

2470 

2570 





925 

.0349 

6.20 

.435 

51 



34.5 

2640 

3395 





1160 



7.00 

.585 

52 



44.5 

2780 

4020 





1370 



7.50 

.685 

53 



49.0 

2860 

4250 





1480 


T 

7 . 70 

.735 

54 



13.4 

2000 

1925 





845 

.0199 

2.50 

.159 

55 



44.2 

2520 

4338 





1785 

.0199 

3.50 

.409 

56 



49.3 

2880 

4213 





1520 

.0349 

7.70 

.715 

57 



49.3 

2890 

4213 


i 


i 

1520 


1 

7.80 

.715 

58 



49.3 

2890 

4213 





1530 


1 

7.80 

i 

.705 
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TABLE II. - Continued. TABULATION OF HEATED AND ISOTHERMAL DATA 


Run 

Mesh 

number 

Voltage, 

V 

Current, 

I, 

amp 

Average 

surface 

tempera- 

ture, 

V 

°r 

Gas 

Inlet 

gas 

tempera- 

ture, 

Tl, 

°R 

Outlet 

gas 

tempera- 

ture, 

°R 

Gas 

flow 

rate, 

V, 

lb/sec 

Inlet 

pressure, 

. p i> 

Ib/sq in. gage 

Total 

static 

pressure 

drop, 

AP t . 

lb/sq in. 

59 

4 

22.4 

1280 

4325 

Nitrogen 

515 

1270 

0.170 





60 



14.8 

1094 

3500 





1070 

.0866 

— 

— 

61 



19.6 

1160 

4200 





1255 

.0866 

— 

— 

62 



13.6 

1010 

3475 





1125 

.0665 

— 

— 

63 



16.6 

1082 

3885 





1255 

.0665 

— 

— 

64 



9.5 

1121 

2365 





785 

.1270 

— 

— 

65 



24.0 

1320 

4485 





1200 

.1270 

— 

— 

66 



26.9 

1396 

4700 





1215 

.1460 

— 

— 

67 



22.8 

1320 

4285 





1115 



— 

— 

68 



16.1 

1280 

3280 





945 



— 

— 

69 



11.8 

1185 

2710 





820 



— 

— 

70 



30.0 

1420 

5065 





1280 



— 

— 

71 



32.0 

1460 

5213 





1335 

T 


— 

— 

72 

3 

16.6 

1000 

3203 

Nitrogen 

515 

865 

0.1460 




73 



23.3 

1080 

2975 





1005 



— 


74 



31.2 

1180 

4725 





1195 



— 


75 



12.6 

960 

2625 





810 

.127 

— 


76 



18.0 

1000 

3420 





945 



— 


77 



23.1 

1044 

4075 





1085 



— 


78 



28.0 

1100 

4588 





1185 



— 


79 



10.9 

920 

2400 





790 

.1070 

— 


80 



17.9 

960 

3515 





990 



— 


81 



24.5 

1040 

4270 





1180 



— 


82 



11.1 

880 

2545 





860 

.0866 

— 


83 



14.2 

880 

3140 





955 



— 


84 



20.0 

920 

4010 





1135 


. 

— 


85 



24.9 

1000 

4485 





1290 


r 

— 


86 



12.4 

840 

2900 





975 

.0665 

— 


87 



16.6 

840 

3700 





1130 





88 



23.0 

940 

4430 




1365 





89 



29.6 

1000 

5173 



1 

1630 


f 

— 


90 

5 

41.0 

3380 

3420 

Hydrogen 

510 

1880 

0.0223 

— 


91 



46.5 

3480 

3700 





2090 



— 


92 



49.3 

3548 

3835 





2223 



— 


93 



37.0 

3140 

3395 





2100 

.0178 

— 


94 



42.4 

3240 ! 

3650 




f 

2257 



— 


95 



46.0 

3340 

3820 



1 

2400 



— 


96 



8.6 

2760 

1085 



535 

685 

.0389 

— 


97 



17.2 

3140 

1765 





885 



— 


98 



26.5 

3320 

2400 





1130 



— 


99 



47.0 

3780 

3488 





1660 



— 


100 



16.8 

3040 

1780 





915 

.0339 

— 


101 



26.5 

3240 

2456 





1200 


! 

— 


102 



32.0 

3380 

2777 





1325 


r 

— 


103 



13.8 

2760 

1630 





865 

.0286 

— 


104 



21.5 

2960 

2225 





1095 

* 

— 


105 

6 

4.5 

3260 

920 

Hydrogen 

560 

645 

0.0389 

— 


106 



7.8 

4040 

1219 





750 



— 


107 



11.1 

4440 

1522 





865 



— 


108 



13.7 

4560 

1780 





955 



— 


109 



28.2 

5060 

3005 



» 

1480 



— 


110 



7.3 

4160 

1130 



540 

740 

.0391 

— 


111 



11.2 

4520 

1510 





865 


t 

— 


112 



19.2 

4800 

2265 





1130 



— 


113 



25.9 

5080 

2790 





1405 



— 


114 



5.7 

3780 

992 




i 

695 

.0339 

— 


115 
1 

1 


11.8 

4320 

1645 





1 

920 

.0339 

— 
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TABLE II. 


Concluded. TABULATION OF HEATED AND ISOTHERMAL DATA 


1 

Run 

Mesh 

Voltage, 

Current, 

Average 

Gas 

Inlet 

Outlet 

Gas 

Inlet 

Total 



number 

V 

I, 

surface 


gas 

gas 

flow 

pressure. 

static 





amp 

tempera- 


tempera- 

tempera- 

rate. 

p i> 

pressure 






ture. 


ture. 

ture. 

W, 

Ib/sq In. gage 

drop. 






T s , 


Tl, 

t 2 . 

lb/sec 

AP t* 

i 





°R 


°R 

°R 



lb/sq In. 


Isothermal 


Hydrogen 520 520 


Hydrogen 520 


Hydrogen 520 


Hydrogen 520 


Nitrogen 520 


Nitrogen 520 520 


Nitrogen 520 | 520 


Nitrogen 520 520 


0.0124 

1.30 

0.0210 

.0178 

2.00 

.0415 

.0232 

2.70 

.0550 

.0285 

3.60 

.0770 

.0338 

4.60 

.0875 

.0391 

5.60 

.1070 

0.0178 

2.00 

0.0915 

.0232 

2.70 

.1450 

.0285 

3.70 

.1970 

.0338 

4.80 

.2570 

.0391 

6.00 

.3170 

0.0178 

2.00 

0.0415 

.0232 

2.70 

.0550 

.0285 

3.50 

.0770 

.0338 

4.60 

.0975 

.0391 

5.80 

.1170 

0.0178 

2.00 

0.0515 

.0232 

2.70 

.0850 

.0285 

3.60 

.1170 

.0338 

4.80 

.1470 

.0391 

5.90 

.1870 

0.0261 

0.70 

0.0080 

.0465 

1.10 

.0170 

.0666 

1.90 

.0310 

.0866 

2.50 

.0500 

.1070 

3.50 

.0670 

.1270 

4.50 

.0870 

.1460 

5.60 

.1040 

0.0465 

1.20 

0.0460 

.0666 

1.90 

.0850 

.0866 

2.70 

.1350 

.1070 

3.60 

.1870 

.1270 

4.65 

.2470 

.1460 

5.90 

.3070 

0.0465 

1.20 

0.0210 

.0866 

2.60 

.0550 

.1070 

3.50 

.0770 

.1270 

4.50 

.0925 

.1460 

5.60 

.1170 

0.0261 

0.70 

0.0090 

.0465 

1.20 

.0230 

.0666 

1.90 

.0410 

.0866 

2.60 

.0700 

.1070 

3.50 

.0990 

.1270 

4.60 

.1280 

.1460 

5.70 

.1620 
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analytical relation that Nusselt number is proportional to the square root of 
the Reynolds number (Nu \/Re (ref. 7)). However, the downstream portion of 
the cylinder where the flow separates defies analysis, and therefore experi- 
mental data must be used to empirically predict the total heat transfer from 
such a surface. The presence of adjacent cylinders affects the boundary layer 
thickness, velocity distribution, and also the nature of the turbulent wake. 

If the cylinders are interwoven into a mesh, the effective heat- transfer area 
is less than the total mesh surface area due to wire overlap. The absence of a 
general correlation for flow through porous media indicates the difficulty in 
defining the geometrical factors that influence both the heat-transfer and 
pressure drop characteristics of the system. 

The usual dimensionless groups were used to correlate the heat- transfer 
data by the relation 


Nu = P[(Re)(Pr)] (13) 

Variable fluid properties and mesh geometry significantly affected the magni- 
tude of these groups. The fluid properties were evaluated at the average bulk, 
film, and surface temperatures to determine which reference temperature best 
represented the variable property effects . The changing geometry that the 
fluid encountered in passing through the mesh was represented by the porosity 
factor included in the definitions of equivalent diameter D e and flow area 
Afp as defined in the section METHOD OF CALCULATION. 

Figures 4(a), (b), and (c) show the correlation of the heat-transfer data 
on a bulk, film, and surface temperature basis, respectively. Evaluation of 
the equilibrium fluid properties at the surface temperature produced the best 
correlation. The physical properties for hydrogen were taken from reference 8 
and for nitrogen from reference 9. The maximum deviation of the data from the 
correlating equation was ±14 percent for the surface temperature correlation 
as compared to ±20 percent and ±29 percent for the film and bulk temperature 
correlations, respectively. 

The equation 

Nu s = 0.462 Re°- 53 Pr°- 40 (14) 

represents the heat-transfer correlation of six helically coiled wire meshes 
for the following range of properties and conditions : 

(1) Porosity of 64 to 72.2 percent 

(2) Wire diameter of 0.020 to 0.035 inch 

(3) Surface temperature of 1400° to 5200° R 

(4) Outlet gas temperature of 600° to 2400° R 

(5) Mass velocity for hydrogen of 0.4 to 3.1 pounds per second per square 

foot and for nitrogen of 4.5 to 10.2 pounds per second per square foot 
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Reynolds number, Re b 

(a) Fluid properties evaluated at bulk temperature. Nu^Pr^ 4 - 0. 303 Re§- 664 



Reynolds number, Re f 

(b) Fluid properties evaluated at film temperature. NUf/Pr^- 4 = 0.50 Re^- 552 . 



(c) Fluid properties evaluated at surface temperature. Nu s /Pr^ 4 = 0.462 Re^ 53 . 


igure 4. - Correlation of he^t-transfer data using equiiiorium fluid properties evaluated at 
bulk, film, and surfdce temperatures. 




(6) Heat flux of 0.5 to 8.3 Btu per second per square inch 

(7) Pressure level of 1 atmosphere 

Practically all of the heat-transfer data presented in the literature for 
a mesh-type configuration were obtained at low-temperature constant property 
conditions. The majority of these data were obtained by the Stanford-ONR 
program on compact heat-transfer surfaces (ref. l). These investigators 
studied several woven wire mesh heat- transfer surfaces of different wire diame- 
ters and porosities. They were able to determine the heat-transfer character- 
istics for each mesh but did not present a general correlation for all of the 
meshes. These data were obtained by a transient test technique whereby the 
mesh test element was heated to a uniform temperature either in a furnace or by 
a stream of hot gas and then immediately subjected to a lower temperature gas 
stream. The time -temperature history of the gas leaving the test element pro- 
vides the information necessary to calculate the heat-transfer characteristics 
of the system. 

In reference 2, the Stanford-ONR data were revised by substituting the 
wire diameter d for the equivalent diameter Dg and applying a correction 
factor for that part of the surface area where the wires overlap and reduce the 
heat-transfer surface. The revised Stanford-ONR data for six meshes resulted 
in one general correlation. 

Reference 3 studied air flow through electrically heated wire tube banks 
for variable property conditions . Their heat-transfer results based on the 
film temperature agreed fairly well with the single wire equation of refer- 
ence 10. 

Figure 5 compares this report's correlation with the correlations of 
references 2, 3, and 10. The Prandtl number was excluded because it was 
similar for all investigations. A direct comparison was not possible because 
of the following differences : 

(1) The fluid velocity for all correlations except reference 10 was based 
on average flow area throughout the mesh as defined by frontal area times 
porosity. For the single wire data of reference 10, the frontal velocity was 
used to characterize the flow. 

(2) References 2, 3, and 10 used wire diameter as the characteristic 
dimension, whereas the equivalent diameter was used for the data of this in- 
vestigation. 

(3) The data of other investigators were obtained at low temperatures 
under relatively constant property conditions as compared to the variable 
property data of this investigation. 

(4) The total surface area was used to evaluate the heat-transfer coef- 
ficient for the data of this investigation because it was not possible to 
accurately define the amount of wire overlap as in reference 2 . 

The good agreement of the single-wire correlation of reference 10 and the 
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Figure 5. - Comparison of other correlations and data with correlation of this investi- 
gation. 


multiwire tube bank data of reference 3 indicates that the effect of adjacent 
wires on heat transfer is not significant if the wires do not engage. The dif- 
ferences in the mesh correlations of reference 2 and this report are indicative 
of the geometry differences between a woven wire mesh and an interwound coil 
mesh. 

The significance of figure 5 lies in the fact that all of the correlations 
have similar slopes approximating 0.5, which indicates that the forced con- 
vective heat transfer for any mesh is primarily a function of the laminar 
boundary layer conditions on the upstream portion of the mesh surface as defined 
analytically by Nu °c C-^/Re (ref. 7). Deviations from this square-root func- 
tion indicate the effect of turbulence created by the different types of mesh 
geometries. Differences in the coefficient C of the various correlating 
equations are due to the geometrical factors such as equivalent diameter, flow 
area, and actual heat- transfer surface area. These factors, which are different 
for each mesh, preclude a universal heat- transfer correlation for all meshes. 

Pressure Drop 

The friction pressure drop, as previously mentioned, is made up of two 
inseparable components (l) a friction component, associated with the surface 
resistance of the wire mesh and (2) an expansion and contraction component, 
associated with the area changes the fluid encounters in passing through the 
mesh. Some investigators (refs. 1, 11, 12, and 13) have used the conventional 
Fanning equation, which is generally used for fully developed flow through 
tubes, to correlate pressure drop data for mesh geometries. This has resulted 
in individual correlations for each mesh. In comparing wire mesh with packed 
beds it became apparent that pressure drop in both is dependent upon the same 
parameters, namely, flow rate, viscosity and density of fluid, porosity, 
packing, size, shape, and surface of the solid. This indicated that a pressure 
drop correlation similar to that used for packed beds might also be applicable 
to wire mesh. Following is a more detailed discussion on the methods used for 
correlating the friction pressure drop data. 
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The conventional Fanning equa- 
tion used to calculate the 
friction-factor data is expressed 
by 


f 

2 


^PgP m D, 


m^e 


(15) 


4LG" 


This equation was used by investi- 
gators in references 1, 11, 12, 
and 13, With the exception of the 
data of reference 3, pressure drop 
data of the other investigations 
has been isothermal. Since this 
report includes pressure drop data 
with and without heat addition for 
hydrogen and without heat addition 
for nitrogen, a variable property 
correlation was used. Figure 6 
shows f/2 plotted .against 
Reynolds number for the heated and 
isothermal data for each of the 
four 3- by 1-inch meshes tested. The mean fluid density was evaluated on a 
bul_k temperature basis (eq. (12)) and the viscosity on a surface temperature 
basis. Other combinations of evaluating fluid density and viscosity for the 
heated data on bulk, film, and surface temperatures were attempted but resulted 
in more scatter. A comparison of the variable property plots of figure 6 with 
similar isothermal plots of references 12 and 13 indicates the following 
similarities exist for meshes of different geometries but within approximately 
the same porosity range: 


Figure 6. - Correlation of average half friction factor with Reynolds 
number using conventional Fanning type correlation. 


(1) The friction factors are of the same order of magnitude 

(2) The slopes of the curves are similar for the same Reynolds number 

(3) A gradual change of slopes occurs between different flow regimes with 

changing Reynolds number 

(4) The curves of the same type mesh sometimes cross 


The conventional Fanning type equation although providing a means of 
evaluating friction factors of individual mesh is limited in use to mesh for 
which pressure drop data has already been experimentally determined. This 
equation provides the only means of comparing the data contained herein with 
that of other investigators. The complexity of the flow through the mesh pre- 
cludes the fact that a general correlation could be obtained by use of the 
conventional Fanning equation. This led to a method of correlation that can be 
easily converted to a modified Fanning type equation, discussed later, and 
which has been successfully used by previous investigators for packed bed 
geometries . 
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Reynolds number, Re s 

Figure 7. - Pressure drop correlation based on method used in reference 14. 



Reference 4 discusses the factors to be considered in correlating parame- 
ters that determine energy losses in packed beds. It further gives a good 
survey of references applicable to packed bed pressure measurements. 


Reference 14 represents the pressure gradient by an equation of the form 


AP 

AL 


F(Re) 

SPDg 


(16) 


where, for laminar flow, AP/AL is proportional to the Reynolds number Re 
and, for turbulent flow, AP/AL is proportional to the square of the Reynolds 
number Re^. In a porous mesh, pressure losses are due to both viscous shear 
and inertia effects associated with laminar and turbulent flow, respectively. 
The pressure gradient according to reference 14 can therefore be represented 
by the sum of the effects of both regimes weighted by the two geometry factors 
a and (3; that is, 

~ = - iL - 3 (a Re + p Re 2 ) (17) 

SPm D e 


Dividing both sides of equation (17) by 


Re |i 2 /gP m D^ 


£PgP m D 2 

ALGp 


a + p Re 


yields 


(18) 


Since a and p are constants dependent upon screen geometry, equation (18) 
is in the form of a straight line (y = mx + b' ) . The advantage of this form of 

friction factor (f = APgp m D 2 /ALGp) is that it is a linear function of the 

Reynolds number and therefore plots as a straight line on an arithmetic scale. 

A plot of the isothermal and heated data as shown in figure 7 indicates a 
straight line can be drawn through the data of each of the four meshes as pre- 
dicted by equation (18). Since the mesh data plot in the order of porosity, 
the addition of a porosity factor to the left side of equation (18) was 
attempted. This should provide a means of plotting the data of all the meshes 
on a single curve. The best factor was found to be a porosity cube term (e 3 ). 
This was determined to be the optimum value, since larger powers of porosity 
although bringing the isothermal data closer, caused the heated data to re- 
spread so it was no longer in the order of porosity. Hence, 


APgp D 2 e 3 

ATif .~ [ = a + p Re s (19) 


where a and jB are 48 and 0.1095, respectively, for the data of this in- 
vestigation. 


Figure 8 shows a plot of APgp m D|e 3 /ALGp s as a function of Reynolds 
number ranging from 200 to 6000. This provides a straight line correlation for 
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Reynolds number, Re $ 


Figure 8. - Pressure drop correlation obtained by modifying method used in reference 14 by porosity 
cube term. 


the data of all the meshes for both heated and isothermal data that has a 
spread of ±25 percent. It should be noted that the viscosity of the fluid is 
evaluated at the mesh surface temperature and the density as defined by equa- 
tion (12) is evaluated at the bulk, temperature. Here again, as with the con- 
ventional Fanning type correlation, the heated data was in best agreement with 
the isothermal data when the fluid properties were evaluated at these temper- 
atures . 

It should be noted that division of both sides of equation (19) by 
Reynolds number transforms this equation into a modified Fanning equation 

dPgpD 

® — = “i. + p = f • (20) 

ALG 2 Re 

where a/Re is the dominating term in the laminar flow regime and (3 is the 
dominating term in the turbulent flow regime and indicates the friction factor 
approaches a constant for turbulent flow. This is clearly shown in figure 6 
where in the turbulent region the friction factors for each mesh approach a 
constant value. 

An equation of the form of equation (20) could have initially been used 
to obtain a correlation; however, a plot of APgp m D e /ZXLG2 against Re for 

each mesh does not plot as a straight line and therefore increases the diffi- 
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Figure 9. - Comparison of present correlating line with data of reference 3. 


culty in obtaining a correlation. Therefore, the approach previously dis- 
cussed was used. 

Figure 9 shows the data of reference 3 for flow of air through banks of 
vertical rods using the correlation of this report (eq. (19)). This indicates 
that a similar correlation might be obtained for different types of mesh. Due 
to geometry differences a single correlation for all mesh would not be ex- 
pected. However, if enough data on different types of mesh were available, 
geometry correction factors could probably be found that would enable a 
universal correlation to be obtained. 


Effect of Screen Bypass 

The correlation of Nusselt number against Reynolds number as discussed in 
the section Heat Transfer has been obtained with the mesh fitting into the 
duct cross-sectional area such that a gap of only 1/64 of an inch (about 
3 percent bypass) exists along each 3-inch side of the mesh. Experimental 
data with hydrogen as the coolant were obtained on one mesh heater of 0.030- 
inch-diameter wire with a bypass area of 25 percent. 

By assuming that all the heat generated in the mesh is transferred to the 


20 



gas passing through the mesh, a heat balance can be made on a mesh with bypass 
as follows : 


1 - W^Pjb^S - Ip) (21) 

The amount of heat transferred to the gas can be represented by the heat- 
transfer equation 


where 


q — hA g ( Tg - Tj-j ) 




( 22 ) 


(23) 


for h, the heat-transfer coefficient, the mesh correlation can be substituted 
as 


h = 0.462 


k s( Ve 

D e\ A fl^s 



(24) 


By combining equations (2l) to (24), the following equation can be written 
for q: 


1 = 


/c p p \ 0,4 k 

0.462 -f- r-2. 

\ k s D e 


f W>e ) 

l A fl^s/ 


0.53 


A b (T h - T-, ) 


1 + 0.462 



0.4 


k srM D e 


,0.53 


(25) 


• D e\ A fl ti s / 


2W M°p,b 


where k s and p s for hydrogen can be represented as approximate functions of 
T g by the following relations: 


= 1.8X10- 7 

rpO .8 

S 

for 500° R > T s 

< 3500° R 

(26) 

= 9 . 1X10“ 8 

tpO . 67 
s 

for 500° R > T 

s 

< 5000° R 

(27) 


For an experimental run the following data were taken W, T s , q, T]_, T 2 , and 
AP^ . Hence all the parameters in equation (25) are known except W^, which 
is then solved for. However, as a means of checking the value of Wjyj- deter- 
mined from equation (25), the total static pressure drop for the mesh was cal- 
culated with equations (ll) and (19), where G is based on the calculated 
value of Wj 4 flowing through the mesh. The value of AP-t calculated from 
equations (ll) and (19) is then compared to the experimental AP-t across the 
mesh and the bypass. The two values of AP-t, measured and calculated, checked 
within ±25 percent, which was within the range of data scatter. For the range 
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of conditions run with a 25 percent bypass 
area, as much as 65 percent of the flow goes 
through the bypass region. 


To illustrate the difficulty encountered 
with bypass, figure 10 compares a set of ex- 
perimental curves for a constant mass flow of 
hydrogen with and without bypass. For greater 
heat-transfer efficiency, a heating element 
should be operated near its maximum tempera- 
ture. From figure 10, it is evident that to 
transfer a given quantity of heat with bypass 
a much higher surface temperature is required 
or for a given surface temperature much less 
heat is transferred. This could cause heater 
burnout if a mesh heater were designed without 
consideration of flow bypass . In order to de- 
sign a heater where the bypass area is no 
longer insignificant (*=3 percent or more), ex- 
perimental bypass flow coefficients are re- 
quired in addition to the heat-transfer and 
pressure drop correlations. 


SUMMARY OF RESULTS 

Heat transfer and friction pressure drop data at 1 atmosphere were ob- 
tained for forced convection of hydrogen and nitrogen through tungsten wire 
mesh electrically heated to surface temperatures up to 5200° R. Outlet gas 
temperatures as high as 2400° R and Reynolds numbers from 200 to 6000 with 
fluid properties evaluated at surface temperatures were obtained. The effect 
of flow bypass, resulting from a mesh heater not filling a flow passage, was 
investigated for a bypass area of 25 percent. The following results were ob- 
tained for the mesh tested in this investigation: 

1. A variable property heat- transfer correlation for helically coiled wire 
mesh operating at surface temperatures up to 5200° R has been obtained based on 
mesh geometry with fluid properties evaluated at the surface temperature. 

2. A pressure drop correlation for both isothermal and heated data has 
been obtained based on mesh geometry with fluid properties evaluated at the 
surface temperature. 

3. Comparison of experimental data with that of other investigators in- 
dicates a universal correlation for all types of wire mesh was not obtained 
because of the difficulty in evaluating the equivalent diameter, surface area, 
and flow area that determine the heat- transfer, pressure drop, and fluid flow 
parameters . 



Figure 10. - Comparison of heat transfer to 
gas with and without bypass from experi- 
mental data. Gas flow rate, 0.0199 pound per 
second; inlet gas temperature, 520° R; gas, 
hydrogen. 
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4. Flow bypassing the mesh heater element can cause serious design 
problems. Experimental data in addition to the heat-transfer and pressure 
drop correlations are required to design a heater with bypass. 


Lewis Research Center , 

National Aeronautics and Space Administration, 
Cleveland, Ohio, April 15, 1965. 
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APPENDIX - SYMBOLS 


A e current flow cross-sectional area, sq ft 
flow area, sq ft 
Af-t frontal area, sq ft 

A s surface area, sq ft 

b length of mesh, ft (except where otherwise noted) 

C constant 

Cp specific heat of gas at constant pressure, Btu/(lb)(°R) 

D mandrel diameter, ft (except where otherwise noted) 

D e equivalent diameter, ft; four times void volume/surface area 
d wire diameter, ft (except where otherwise noted) 

F function 

f average conventional friction factor 

f ' average modified friction factor 

G mass velocity through mesh, lb/(sec)(sq ft) 

g acceleration due to gravity, 32.2 ft/sec^ 

h average heat-transfer coefficient, Btu/(sec)(sq ft)(°R) 

I current, amps 

k thermal conductivity of gas, Btu/(ft)(sec)(°R) 

L thickness of mesh in direction of flow (characteristic length), ft 

N number of parallel coils 

Nu Nusselt number, hD g /k 

P pressure, lb/sq ft 

Pp inlet static pressure, lb/sq ft 

P2 outlet static pressure, lb/sq ft 

AP friction pressure drop, lb/ sq ft 
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AP-t total static pressure drop, lb/ sq ft 
Pr Prandtl number, Cpp/k 

p coil pitch, ft/turn 

q. rate of heat transfer to gas, Btu/sec 

R gas constant, (ft)(lb)/(lb)(°R) 

Re Reynolds number, D e G/p 

S total wire length of helical coil, ft 

T temperature, °R 

T-j-, average bulk temperature, (T^ + Tg)/2, °R 

average bulk temperature of unmixed gas passing through screen, 
T-l + Tg/2, °R 

T s average surface temperature, °R 

Tp inlet gas temperature, °R 

Tg outlet gas temperature, °R 

Tg outlet temperature of unmixed gas passing through screen, °R 

V voltage 

W gas flow rate, lb/sec 

Wjfl gas flow rate through mesh, lb/sec 

a experimental screen geometry factor (viscous flow) 

3 experimental screen geometry factor (turbulent flow) 

e screen porosity 

£ resistivity of tungsten, ohm- ft 

H absolute viscosity of gas, lb/(sec)(ft) 

mean gas density, (Pp + Pg)/R(Tp + Tg), lb/cu ft 
inlet gas density, lb/cu ft 
P 2 outlet gas density, lb/cu ft 
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Subscripts: 
b bulk 

f film 

s surface 
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